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An analysis of the optimal control by blowing and suction in order to generate stream- 
wise velocity streaks is presented. The problem is examined using an iterative process 
that employs the Parabolized Stability Equations for an incompressible fluid along with 
its adjoint equations. In particular, distributions of blowing and suction are computed for 
both the normal and tangential velocity perturbations for various choices of parameters. 


Nomenclature 


E energy norm 

G energy growth ratio 

L re f reference length 

p pressure disturbance 

U mean-flow streamwise velocity 
U re f reference velocity 
u streamwise velocity disturbance 
V mean-flow normal velocity 

V w normal velocity of blowing or suction at the wall 
v normal velocity disturbance 

W w tangential velocity of blowing or suction at the wall 

w spanwise velocity disturbance 

/ 3 spanwise wave number 

e \/n/U re fL re f 

v kinematic viscosity 

p density 

Superscript 
T transpose 
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Subscript 

exit exit conditions 
in start point 
out end point 
ref reference 
x d/dx 

V d/dy 


I. Introduction 

The performance of low-pressure turbines (LPT’s) is strongly affected by flow separation. There is 
a possibility of delaying the boundary layer separation by tripping the boundary layer with the use of 
generators placed on the surface of the airfoil. Usually, a trial-and-error method is used to determine an 
appropriate placement of the control elements. However, this approach is time consuming and expensive. 

One of the possible scenarios of early transition is associated with the development of velocity streaks. 
Due to a secondary instability mechanism (associated with the presence of free stream turbulence), velocity 
streaks in the boundary layer flow may lead to bypass transition at low Reynolds numbers. This earlier 
transition to turbulence can prevent the flow separation. To generate the streaks, one may introduce a 
system of counter-rotating streamwise vortices that decay downstream and result in a streaky structure in 
the boundary layer. 

Analysis of optimal streamwise vortices within the scope of the transient growth theory 1 indicated that 
a favorable pressure gradient suppresses the transient growth mechanism. Therefore, a row of localized 
generators placed on the wall will be unlikely to generate the optimal perturbation. However, there is also 
a possibility that the same streaky structures can be generated with the help of distributed generators. 

In Ref. 2, a model of optimal control of the algebraically growing perturbations by blowing and suction at 
the wall was proposed. The model provided a tool to find the distribution of blowing and suction at the wall 
that suppresses the streaky structures. Because the theoretical model 2 is based on linearized equations, it 
can also be used to optimize the distributed generators in order to generate the streaky structures equivalent 
to those generated by the optimal perturbation. 

In the present work, the distributed generators are emulated by the normal and tangential velocity 
perturbations at the wall. A parametric study of these generators has been performed for a flat plate 
without a pressure gradient and for a flat plate with a favorable pressure gradient. 

II. Governing Equations 

We consider steady three-dimensional disturbances in an incompressible two-dimensional boundary layer. 
The streamwise coordinate, x, is made nondimensional using a reference length scale L re f, while the 
normal and spanwise coordinates, y and z respectively, are made nondimensional using the length scale 
\J vL re f / U re f , where v is viscosity and U re f is a reference velocity. With the help of a small parameter e 
defined as e = i , /U re fL re f, the following scalings are assumed for the streamwise, normal and spanwise 
velocity disturbances, u, v and w respectively, and the pressure disturbance, p: 

u ~ U re f, V^eUref, W~eU re f, p ~ £ 2 pU re f (1) 

This choice of scaling of the linearized Navier-Stokes equations along with the neglect of curvature effects leads 
to the governing equations for Gortler instability with Gortler number G = 0. The solution is assumed to be 
periodic in the spanwise direction z. It is therefore possible to assume solutions of the form u = u(x, y) cos (3 z, 
v = v(x, y) cos j3z, w = w(x, y) sin/fe and p = p(x, y) cos /3z where j3 is the spanwise wave number. The 
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governing equations are given in dimensionless form as: 


u x + v y + /3w = 0 

(2) 

( Uu) X + Vu y + VUy = Uyy ~ /U M 

(3) 

(' UV + VU)x + (2 Vv)y + /3VW + Py = Vyy — fl 2 V 

(4) 

( UW) X + ( Vw)y - PP = Wyy ~ ffP W 

(5) 


where U(x,y) is the mean-flow streamwise velocity and V(x,y) is the mean-flow normal velocity. U is scaled 
with U re f and V is scaled with eU re f. 

Since we investigate two types of distributed generators, there are two sets of boundary conditions. The 
first set of boundary conditions is associated with the normal velocity perturbation and is given as: 

y = 0 : u = w = 0, v = V w and y — > oo : u, w, p = 0 (6) 

where V w is the normal velocity of blowing or suction at the wall (y = 0). The second set of boundary 
conditions is associated with the tangential velocity perturbation and is given as: 

y = 0: u = v = 0, w = W w and y — » oo : u,w,p = 0 (7) 

where W w is the tangential velocity of blowing or suction at the wall. The system of equations (2-5) can 
be solved subject to the boundary conditions (6) or (7) with a prescribed initial velocity perturbation at 
x = Xq. 

Introducing the vector f = (u,v,w,p) T , where the superscript “T” means transpose, the governing 
equations (2-5) can be recast as 

(-A-f)x = Bof + Bify + (8) 

where A, Bo, Bi, and B 2 are 4x4 matrices given in the appendix as well as in Ref. 1 and Ref. 3. 

III. Optimal Control: Model 


A. Model 

We are interested in computing the distribution of blowing and suction that can generate a streamwise 
velocity streak equivalent to that which will be generated by an initial optimal perturbation. Since the 
problem is linear, this distribution can be determined using the optimal control model of Cathalifaud and 
Luchini. 2 

Given an optimal initial disturbance at Xq, the method determines a control by blowing and suction 
that suppresses the streamwise streaks. The method of optimal control is based on the minimization of 
an objective functional subject to constraints. Using standard mathematical techniques from control the- 
ory, 3 4 a system of differential equations (direct and adjoint problems) along with conditions of optimality is 
formulated. The solution gives the optimal control and is computed using an iterative procedure. 

For the sake of clarity, we recapitulate the formulation of the optimal control model. In this problem, 
f is the vector containing the state variables, V w (or W w ) is the control variable and the constraints are 
given by the direct problem (Eq. 8) along with its boundary conditions and initial condition. Additionally, 
£= ( a,b,c,d) T is the vector containing the Lagrange multipliers and £ is the Lagrange multiplier associated 
with the wall receptivity. Although we continue the model description using the control variable V w , the 
model is valid for W w as well. The objective functional to be minimized is given as 

1 f°° 

J(f,V w ) = -j u 2 (x out ,y) dy (9) 
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Minimization of this functional is equivalent to minimizing the energy of the disturbance at the output 
cross-section. Note that this objective functional does not take the “cost” of the control into account. 

The Lagrangian functional is defined as 


r%out rOO pXout 

£(f,V w ,€,0 = J(f,V w )- / £-[(Af) x -B 0 f — Bify — B 2 f TO ] dy dx- / C{x)-{v{x,o)-V w {xj) dx 

” Win ^ 0 "Win 

( 1 °) 

and the equations of the adjoint problem can be found from the condition of extremality 


sc = Um £(f + ^f,K,,g,Q-£(f,y m ,$,C) 

r — *0 r 

where <5f is a variation of f. The equations for the adjoint problem are 


( 11 ) 


a X + Ub X + V C X + (Vb)y + byy ~ ^ 6 — 0 

U C X T 2 PCy “t“ Cyy (3 C Uyb dy = 0 

ud x + V d y + d yy — ft 2 d — (3{a + V c) =0 


The boundary conditions are 


c y + /3d = 0 


(12) 

(13) 

(14) 

(15) 


y = 0 : a = b = c = 0 and y — > oo : a + c y + 2 cV = 0 (16) 

and the initial condition of the adjoint equations, u a (x ou t,y ) is given by the solution to the direct problem 
at the end point, u{x ou t,y). Furthermore, the optimality condition is given by 

y = 0 : a + c y = 0 (17) 


B. Iterative Procedure 

An iterative procedure is used to numerically determine the optimal control. The algorithm is given as: 

1. The forward problem using the direct equations is solved from x = Xi n to x = x ou t subject to the 
boundary conditions and the initial condition (the optimal perturbation) . 

2. The objective functional is computed at the end of the forward iteration, and it is compared to the 
value of the objective functional found at the end of the first iteration. If the ratio is less than a prescribed 
value, the iterations stop. Otherwise, the initial condition for the backward problem is found from the solu- 
tion to the forward problem at x = x ou t- 

3. The backward problem using the adjoint equations is solved from x = x ou t to x = Xi n subject to the 
boundary conditions and the initial condition for the adjoint problem. 

4. A new initial condition and a new boundary condition for V w is computed, and the forward problem 
is iterated. The loop continues until the optimization is considered to have converged. 

Given an initial optimal perturbation, this iterative scheme determines the optimal control by blowing and 
suction to suppress the streamwise velocity streak that would have been generated by the initial disturbance 
(if there was no control). Again, since the problem is linear, this same optimal control by blowing and 
suction can be used to generate this same velocity streak in the case of zero initial disturbance. 

4 of 13 


American Institute of Aeronautics and Astronautics 



IV. Results 


A. Flat Plate Without a Pressure Gradient 

We consider a flat plate of length L in a flow having free stream velocity U^. Therefore, this case has 
the scalings L re f = L and U re f = U a 0 . The flow control area is distributed over the interval [xo, 1], where 
Xo, the location of the optimal inflow disturbance, is varied. For each pair of the parameters (3 and xo, a 
transiently growing optimal inflow was determined using the model described by Andersson, et.al., 3 Luchini 5 
and Tumin, et.al. 1 

The optimal inflow disturbance is found in terms of the energy growth ratio G = E out /Ei n , where Ei n is 
an input energy norm and E out is an output energy norm. Employing the fact that the optimal disturbances 
at the input are represented by streamwise vortices with corresponding streamwise streaks at the output, we 
use the following energy norms of Ref. 1 and Ref. 5: 

rUmax rUmax 

E in = e 2 / ( v 2 + w 2 ) dy and E out = / u 2 dy (18) 

Jo Jo 

As an example, the streamwise velocity perturbation at x = 1.0 is shown for (3 = 0.4 and Xq = 0.25 in 
Fig. 1. Figure 2 shows the optimal velocity perturbations at xo = 0.25 that corresponds to the streamwise 
velocity perturbation at x = 1.0 shown in Fig. 1. 




Figure 1. Streamwise velocity perturbation at x = 1.0 Figure 2. Optimal velocity perturbations for /3 = 0.4 
for 3 = 0.4 and xo = 0.25. at xo = 0.25 corresponding to the streamwise velocity 

perturbation at x = 1.0 shown in Fig. 1. 

After the optimal inflow was computed, the iterative procedure described in Section III was applied to 
calculate the wall velocity distribution that would suppress the velocity streaks at x = 1.0. Due to the linear 
character of the problem, this velocity distribution is also the optimal distribution for the generation of a 
streamwise velocity streak that has the same properties as in the case of an optimal inflow perturbation. 

Figures 3-6 depict the blowing and suction distribution of V w , the normal velocity at the wall, that will 
lead to optimal streaky structures for various combinations of the parameters (3 and Xq. Figure 3 shows V w 
for xo = 0.1 for four choices of [3 {(3 = 0.2, 0.4, 0.6 and 0.8). For these four values of (3, Fig. 4 shows V w for 
xo = 0.25, Fig. 5 shows V w for xo = 0.4 and Fig. 6 shows V w for xo = 0.55. 

In each figure, one can see that for optimal control, a large region of blowing exists at the start of the 
flow control interval. For each choice of Xo, the magnitude of V w at Xo increases as f3 increases. Also, for 
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each choice of (3, the magnitude of V w at Xq increases as Xo increases. After the interval of blowing, there is 
a region of suction. The location of this blowing-suction transition within the flow control interval remains 
roughly the same for a choice of (3. Within this region of suction, the maximum magnitude of V w is associated 
with the largest value of (3 for a given xq', also, this magnitude increases as Xq increases. Beyond this region 
of suction there are alternating smaller regions of blowing and suction. The number of these regions increases 
as xo increases. Furthermore, because of the scalings, the physical size of V w is very small. 

Figures 7-10 depict the blowing and suction distribution of W w , the tangential velocity at the wall, that 
will lead to optimal streaky structures for various combinations of the parameters f3 and Xo ■ Figure 7 shows 
W w for Xo = 0.1 for four choices of (3 ((3 = 0.2, 0.4, 0.6 and 0.8). For these four values of (3 , Fig. 8 shows 
W w for xq = 0.25, Fig. 9 shows W w for xq = 0.4 and Fig. 10 shows W w for xq = 0.55. 

Statements similar to those made about Figs. 3-6 can be written regarding the relative magnitudes of 
W w for fixed values of Xo and /3. One major difference is that the first region in the flow control interval 
is one of suction. This region is followed by a region of blowing, and then there are alternating regions of 
suction and blowing, with the number of regions increasing as xo increases. As with V Wl because of the 
scalings, the physical size of W w is very small. 

B. Flat Plate With a Pressure Gradient 

Optimal perturbations in a boundary layer with a pressure gradient were considered in Ref. 1. The pressure 
distribution corresponded to the model designed at the NASA Glenn Research Center that was used to study 
the flow over a LPT blade. This case has the scalings L re f = L s , where L s is the suction surface length 
and U re f = U exit , where U ex u is the nominal exit free-stream velocity based on the inviscid solution. 1,6 The 
optimal inflow disturbance was applied at Xo/L s = 0.125 and Xq / L s = 0.25 with the optimal energy growth 
occurring at x/L s = 0.5. For each choice of Xo/L s , four values of (3 were considered (/ 3 = 0.3, 0.5, 0.7 and 
0.9). 

We again apply the optimization technique to find the blowing and suction distributions for both the 
normal and tangential velocities at the wall in order to generate the same streaky structures as in the case 
of optimal perturbations. Figures 11 and 12 show V w for several choices of (3 where Xo/L a = 0.125 and 0.25 
respectively. For these same values of (3, Figs. 13 and 14 show W w for Xo/L s = 0.125 and 0.25 respectively. 

The general behavior of the relative magnitudes of V w and W w for choices of Xq and (3 remain consistent 
with that shown in Figs. 3-6 and Figs. 7 - 10. However, in the case of a flat plate with a pressure gradient, 
the first region of V w in the flow control interval is a region of suction. This region is followed by a region of 
blowing, and again these regions are followed by alternating regions of suction and blowing, with the number 
of regions increasing as Xo increases. Additionally, for the case of a flat plate with a pressure gradient, the 
first region of W w in the flow control interval is one of blowing. As with V w , this region is followed by 
alternating regions of suction and blowing. 


V. Conclusions 

In this paper we used an iterative method to compute the optimal control by blowing and suction at the 
wall that minimizes the energy of perturbations when the initial perturbation is optimal. Due to the linear 
nature of the governing equations, these same distributions can be used to generate the streamwise velocity 
streaks that would have been generated from the optimal disturbance. The magnitude of the normal and 
tangential velocity perturbations that can generate the velocity streaks is very small. Therefore, unlike the 
case of localized generators, distributed generators may be a good candidate to generate the velocity streaks 
that lead to bypass transition. 
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Appendix 


The matrices of Eq. (8) are given as follows: 
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Figure 3. The blowing and suction distribution for the normal velocity at the wall for the flow control interval 
[0.1, 1] for four choices of / 3 . 



Figure 4. The blowing and suction distribution for the normal velocity at the wall for the flow control interval 
[0.25, 1] for four choices of ( 3 . 
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Figure 5. The blowing and suction distribution for the normal velocity at the wall for the flow control interval 
[0.4, 1] for four choices of ( 3 . 



Figure 6. The blowing and suction distribution for the normal velocity at the wall for the flow control interval 
[0.55, 1] for four choices of ( 3 . 
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Figure 7. The blowing and suction distribution for the tangential velocity at the wall for the flow control 
interval [0.1, 1] for four choices of ( 3 . 



Figure 8. The blowing and suction distribution for the tangential velocity at the wall for the flow control 
interval [0.25, 1] for four choices of ( 3 . 
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X 


Figure 9. The blowing and suction distribution for the tangential velocity at the wall for the flow control 
interval [0.4, 1] for four choices of ( 3 . 



Figure 10. The blowing and suction distribution for the tangential velocity at the wall for the flow control 
interval [0.55, 1] for four choices of ( 3 . 
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Figure 11. The blowing and suction distribution for the normal velocity at the wall for the flow control interval 
[0.125,0.5] for four choices of ( 3 . 



Figure 12. The blowing and suction distribution for the normal velocity at the wall for the flow control interval 
[0.25, 0.5] for four choices of ( 3 . 
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Figure 13. The blowing and suction distribution for the tangential velocity at the wall for the flow control 
interval [0.125,0.5] for four choices of ( 3 . 



Figure 14. The blowing and suction distribution for the tangential velocity at the wall for the flow control 
interval [0.25, 0.5] for four choices of / 3 . 
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